Crossing the Phantom Divide Barrier with Scalar Tensor Theories 
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There is accumulating observational evidence (based on Snla data) that the dark energy equation 
of state parameter w may be evolving with time and crossing the phantom divide barrier w = 
— 1 at recent times. The confirmation of these indications by future data would indicate that 
minimally coupled quintessence can not reproduce the observed expansion rate H(z) for any scalar 
field potential. Here we explicitly demonstrate that scalar tensor theories of gravity (extended 
quintessence) can predict crossing of the phantom divide barrier. We reconstruct phenomenologically 
viable scalar-tensor potentials F(<&) and t/($) that can reproduce a polynomial best fit expansion 
rate H(z) and the corresponding dark energy equation of state parameter w(z) crossing the w = —1 
line. The form of the reconstructed scalar tensor potentials is severely constrained but is not uniquely 
determined. This is due to observational uncertainties in the form of H(z) and also because a single 
observed function H(z) does not suffice for the reconstruction of the two potential functions F(<i>) 
and [/($). 



PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 



parameter H(z) as a function of the redshift z given by 
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Detailed observations of the relation between luminoc- 
ity distance and redshift for extragalactic Type la Super- 
novae (Snla) indicatesplQ that the universe has entered 
a phase of accelerating expansion (the scale factor obeys 
a > 0). In addition, a diverse set of other cosmologi- 
cal data which includes large scale redshift surveys Q 
and measurements of the cosmic microwave background 
(CMB) temperature fluctuations spectrum £| has shown 
that the spatial geometry of the universe is flat but there 
is not enough matter density to justify this flatness. Thus 
there is an additional cosmological component required 
to justify the observed flatness. This component should 
have repulsive gravitational properties to justify the ob- 
served present accelerated expansion. Such properties 
can either be due to a modified theory of gravity 5J or, in 
the context of standard general relativity, to the existence 
of a smooth energy component with negative pressure 
termed 'dark energy' 0, Q- This component is usually 
described by an equation of state parameter w = - (the 
ratio of the homogeneous dark energy pressure p over 
the energy density p). For cosmic acceleration, a value 
of w < — | is required as indicated by the Friedmann 
equation 
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(p + 3p) 



(1.1) 



The cosmological constant (w = —1) is the simplest 
viable candidate for dark energy. It predicts an expansion 
history of the universe which is described by a Hubble 



= ij 2 [n 0m (i + zf + (i - n 0m )] 

(1.2) 

where flatness has been imposed and Qom = Po/pc is 
the single free parameter of this simplest data consistent 
parametrization (LCDM). Such a model has two impor- 
tant disadvantages: 

• It requires extreme fine tuning of the value of the 
cosmological constant A for the accelerating expan- 
sion to start at around the present cosmological 
time (coincidence problem). 

• It does not provide the best possible fit to recent 
Snla data. 

In particular, recent analyses 0, |j| 0, El of the lat- 
est and most reliable Snla dataset (the Gold datasetQ) 
have indicated that significantly better fits compared to 
LCDM are obtained by allowing for a redshift depen- 
dent equation of state parameter. Extensive analyses 
of such parametrizations of H(z) have shown that the 
parametrizations that allow crossing of the w = —1 line 
(known also as Phantom Divide Line - PDL) provide sig- 
nificantly better fits to the data. It is therefore important 
to construct physical models that allow for a rcdshift- 
dependent w that crosses the PDL. It has been shown 
however that this is a highly non- trivial task[H[ll[l|. 

The simplest approach towards constructing a physi- 
cal model for dark energy is to associate it with a homo- 
geneous minimally coupled scalar field (j> with negative 
pressure whose dynamics is determined by a potential 
properly chosen so that the energy density of <f> comes to 
dominate the universe at the present time. Such models 
are described by Lagrangians of the form 
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where the upper (lower) sign corresponds to a 
quintessence [l5| (phantom !]!!) fi eld. It should be noted 
however that phantom fields with negative kinetic 
term (w < —1), violate the strong energy condition, 
the null energy condition (p + p > 0) and maybe phys- 
ically unstable. They also lead to a future 'Big Rip' 
singularity where all bound systems get dissociated due 
to the increasing repulsive gravitational force of phan- 
tom energy^?]- The phantom instability however can be 
cured in extended gravity theories [l^. The equation of 
state parameter corresponding to is 



rate H(z) that violates the inequality 



w = - = 



±\fi + V{4>) 



(1.4) 



For quintessence (phantom) models with V{<j>) > 
(V{<f) < 0) the parameter w remains in the range 
— 1 < w < 1. For an arbitrary sign of V(<p) the 
above restriction does not apply but it is still impos- 
sible for w to cross the PDL w = —1 in a continous 
manner. The reason is that for w = — 1 a zero ki- 
netic term ±</> 2 is required and the continous transition 
from w < — 1 to to > - 1 (or vice versa) would require 
a change of sign of the kinetic term. The sign of this 
term however is fixed in both quintessence and phantom 
models. This difficulty in crossing the PDL w = — 1 
could play an important role in identifying the correct 
model for dark energy in view of the fact that data favor 
w ~ — 1 and furthermore parametrizations of w(z) where 
the PDL is crossed appear to be favored over the cosmo- 
logical constant w = — 1. Even for generalized k-essence 
Lagrangians ^1[2(J of a minimally coupled scalar field eg 



V{4>) 



(1.5) 



it has been shown [12] to be impossible to obtain cross- 
ing of the PDL. Multiple field Lagrangians (combinations 
of phantom with quintessence fields [2lLl22l 123 . |24 | ) have 
been shown to in principle achieve PDL crossing but such 
models are complicated and without clear physical mo- 
tivation (but see j2|| for an interesting physically moti- 
vated model). 

The obvious class of theories that could lead to a solu- 
tion of the above described problem is the non-minimally 
coupled scalar fields. Such theories are realized in a uni- 
verse where gravity is described by a scalar-tensor theory 
and their study is well motivated for two reasons: 

1. A scalar-tensor theory of gravity is predicted by all 
fundamental quantum theories that involve extra 
dimensions. Such are all known theories that at- 
tempt to unify gravity with the other interactions 
(eg supergravity (SUGRA), M-theory etc). 

2. As shown in the following sections, quintessence 
scalar fields emerging from scalar tensor theories 
(extended quintessence) can predict an expansion 



d(H(zf/Hl) 
dz 



(1.6) 



It is easy to show that the inequality i|l.(jfl is equiv- 
alent to 



w(z) 
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(seeegRef. pU for a derivation of w(z) in the form 
of equation (|1.7|) ). Thus, violation of the inequality 
11.61) is equivalent to crossing the PDL w = — 1 and 
is favored by the Gold Snla dataset [HIIIEIEI 
l29l l3fj| . The inequality (|1.6[) can not be violated 
in minimally coupled quintessence theories (see eq. 
(|2~T2l ). 

The usual approachpj |M El HI S 13 m comparing 
quintessence models with observations is to start with a 
well defined theoretical model Lagrangian, identify the 
predicted form of H (z) and compare with observational 
data to identify the consistency and the quality of fit of 
the model. 

This approach is not particularly efficient in view of 
the infinite number of possible model Lagrangians that 
may be considered. An alternative, more efficient ap- 
proach is to start from the best fit parametrization H{z) 
obtained directly from data and use this H(z) to recon- 
struct the corresponding theoretical model Lagrang ian. 
This later approach was pioneered in Refs |35l l36t l37j| 
and has been further developed for the cases of both min- 
imally coupled quintessence [H EH E3 and scalar tensor 
theories |4ll E^ (extended quintessence [H El E l E ^ L 
Extensions of this approach have recently also been ap- 
plied to f{R) generalized gravity theoriespTj. However, 
despite the high quality of the data of the Gold dataset 
that allows a fairly reliable determination of H(z) (espe- 
cially at redshifts up to z ~ 1) no attempt has been made 
to reconstruct quintessence or extended quintessence La- 
grangians from the H(z) that best fits the Gold dataset. 
This task is undertaken in the present study. 

We consider a simple polynomial parametrization of 
H(z ) a nd fit it to the Gold dataset following Refs. 
Ill El El- The best fit form of H{z) is found to vi- 
olate the inequality (|1.6f> ie to cross the PDL. Thus it is 
not consistent with minimally coupled quintessence. We 
thus use it to reconstruct a class of extended quintessence 
scalar-tensor Lagrangian which involves two scalar func- 
tions F(<5) and t/($). 

The structure of the paper is the following: In section 
II we present in some detail the reconstruction method 
described above for the case of extended quintessence and 
derive the main equations that relate the scalar-tensor 
potentials F(&) and U($>) with the observed Hubble pa- 
rameter H(z) and other experimental constraints. In sec- 
tion III we apply this technique to reconstruct a class 
of F(Q) and {/($) from the actual best fit polynomial 
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H(z) obtained from the Gold dataset. The reconstructed 
potentials are consistent with other observational con- 
straints (solar system tests and Cavendish type experi- 
ments). We show that in contrast to minimally coupled 
quintessence, this reconstruction is possible for extended 
quintessence. However, even though the reconstructed 
potentials F(<&) and [/(<&) are severely constrained by 
the form of H(z) and other consistency requirements, 
they are not uniquely determined. Finally, in section 
IV we review the prospects for a more constrained recon- 
struction of scalar-tensor theories which may significantly 
contribute towards the identification of the fundamental 
theory realized in Nature. 

II. TECHNIQUE FOR RECONSTRUCTING 
THE SCALAR-TENSOR LAGRANGIAN 

The scalar-tensor Lagrangian is a generalization of the 
general relativistic Lagrangian of the form 

F($) Z($) 

C= -Li R Li g^d^dv* -[/($) + £mbPm;9^] ■ 

(2.1) 

where we have set 8irG = 1 (Fq = 1) and C m repre- 
sents the matter fields and does not depend on <E> so that 
the weak equivalence principle is satisfied. One of the 
degrees of freedom F($) and Z(&) can be absorbed by 
a redefinition of $. The following parametrizations are 
commonly used 

• F($) -> $, Z($) -> ^r 1 is the Brans-Dickc (BD) 

parametrization where lo(<&) = rajyg^p is the BD 
parameter. 

• Z($) — > ±1. The case Z — » — 1 corresponds plj 
to a BD parameter —3/2 < ojq < which contra- 
dicts solar system tests [48] for any {/($) allowing 
cosmological evolution of <I>. 

• 9w - 9% = 9v> and $ -> <p : (ff)* = 



3 / d ln ^W ^ _|_ 2F(2) ■ ^ ms transformation corre- 
sponds to the Einstein frame as opposed to the 
original Jordan frame of equation (|2.1|) . In the 
Einstein frame the kinetic terms of the graviton 
and the scalar field are diagonalized and the math- 
ematical analysis of the theory is simplified at the 
expence of introducing an explicit coupling of the 
scalar field with matter. An advantage of the origi- 
nal Jordan frame is that the various physical quan- 
tities are those measured in experiments. 

In what follows we choose the Jordan frame and the 
parametrization Z — ► 1. Under this assumption, F(&) 
should satisfy the following constraints 

• F($>) > so that gravitons carry positive 
energy 41] . 

• Uq 1 = (dF/dQ) 2 , < 4 x 10~ 4 from solar sys- 
tem measurements |48j (the subscript o denotes the 
present time). 

Assuming a homogeneous $ and varying the action cor- 
responding to H2.1fl in background of a flat FRW metric 

ds 2 = -dt 2 + a 2 (t)(dr 2 + r 2 (d9 2 + sin 2 6 d(f> 2 )) (2.2) 
we find the coupled system of equations [42| 

3F-H 2 = p+~$ 2 -3H -F + U (2.3) 
-2FH = (p + p) + $ 2 + F - H ■ F (2.4) 

where we have assumed the presence of a perfect fluid 
(p,p). Eliminating $ 2 from 1)2.4(1 . setting 

q(z) ee H(z) 2 /H 2 (2.5) 

and rescaling U — > U ■ Hq while expressing in terms of 
redshift z we obtain 
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q' 4 
~2q 1 + z 
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1 + z (1 + z) 2 (l + z) 2 q 2 
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r 



(2.6) 
(2.7) 



where the prime ' denotes differentiation with respect to 
redshift (4^) and we have assigned properties of matter 
Op = 0, Oom = -fjr 5 -) to the perfect fluid. 

Given the form of H (z) from observations, equations 
(|2.ti|) and l|2.7|l may be used with boundary conditions 
F(0) = 1, F'(0) = and $(0) = to reconstruct sets of 



F(z), U(z), $(z) which predict the given form of H(z) 
and are consistent with solar system tests. The system 
((2.6(1 and 1(2. 7JI may take a more convenient form by set- 
ting 

m - (2.8) 
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It then takes the form 



U(z) = + \</(l + z) 4 f-~(l + z) a n 0m (2.9) 

$'(z) 2 = -6/-6(l + z)/'-f (l + z )2 / '-(l + z) 2 /"-3fi 0m ^ ±£ (2.10) 

I 



with boundary conditions $(0) = 0, /(0) = 1 and /'(0) = 
-2. 

The case of minimally coupled $ is regained for F(z) = 
1 (/ = 1/(1 + z) 2 ). Indeed, setting f(z) = 1/(1 + z) 2 in 
Il2^1l- (I2~T71 we find 

tf(z) = 3<7-|(l + z)-^(l + z) 3 r! 0m (2.11) 
4>' 2 (z) = * J q' - 3r> 0m (l + zf) (2.12) 

which are in agreement with previous studies [35l l38j] re- 
constructing minimally coupled Lagrangians. 

The reconstruction of [/($) may proceed in the mini- 
mally coupled case by finding §>(z) from 12.12f) (setting eg 
$(0) = 0), inverting for z($) and substituting in l|2.11|l to 
find U(&). The inequality IfLGfi (preventing PDL cross- 
ing), valid for the minimally coupled case is recovered 
from (f^TTSf) since <I>' 2 > 0. 

As a simple application we may consider the H(z) in- 
duced by a cosmological constant which implies 

q(z) = Q A + n am (l + z) 3 (2.13) 
In this case H2.llfH2.12fl reduce to 

U(z) = m A = PA (2.14) 
$'(z) = (2.15) 

as expected. Thus, the reconstruction leads to a uniquely 
defined potential in the minimally coupled case if q' > 

3ft m(l + z) 2 . 

If on the other hand there is a redshift range where 
q' < 3fl 0m (l + z) 2 , an F(z) + \{f ± 1/(1 + zf) is re- 
quired to keep $' 2 > in (|2TTU|) . Since q' < 3Q 0m (l + z) 2 
implies superacceleration which can not be supported 
merely by the repulsive gravity of the potential U , a mod- 
ified strength of gravity will be required (F(z) < 1). Such 
an F(z) is also constrained to obey F(z) > for all red- 
shifts and F'(z = 0) ~ 0. It is therefore possible that 
no such F(z) exists so that for a given q{z) the three 
requirements 

$'(z) 2 > (2.16) 
F(z) > (2.17) 
F'{z = 0) = (2.18) 



are fulfilled. In fact, for q{z) corresponding to a cosmo- 
logical constant (equation (|2.13fl 'l and U(z) = it has 
been shown that no acceptable function F(z) exists |42^ 
because F(z) becomes negative at relatively low z. If 
an acceptable F(z) is found to exist then it will not be 
uniquely defined because a perturbed F{z) will also lead 
to a positive <I>' 2 leading to another acceptable solution. 

The best case scenario for the existence of an accept- 
able F(z) is the case when $'(z) = ie when F(z) does 
not have to balance the attractive gravity of <& /2 . If no 
acceptable F(z) exists for $' = (ie if F(z) < for 
z > z c ) then there will be no acceptable F(z) for any 
<f>' 2 > (F(z) will get even more negative to balance the 
larger attractive gravity of $ /2 ). Thus, in order to see if 
an acceptable scalar tensor theory exists for a given q(z) 
we simply have to set $' = in equation (|2.1Uf) . solve 
for f(z) with /(0) = 1, f(0) = -2 and see if the solu- 
tion obeys f(z) > 0. If f(z) < for some redshift range 
then it will remain so for any $' 2 > and therefore no 
acceptable scalar tensor Lagrangian exists predicting the 
given q(z). If on the other hand f(z) > for all red- 
shifts of interest, the solution is acceptable and we may 
start increasing <i>' 2 (trying different functional forms) 
until f(z) < for some redshift range. Thus we may 
identify a set of observationally acceptable scalar tensor 
Lagrangians that predict the given form of q(z). This 
procedure will be applied in the next section to iden- 
tify scalar-tensor Lagrangians that predict the particular 
form of q(z) that consists one of the best fits to the Gold 
dataset. 



III. FROM THE GOLD DATASET TO THE 
FUNDAMENTAL THEORY 

A wide range of expansion history parametrizations 
H(z) with a flat priori [HJ[II1 has been fitted to the 
Gold and earlier Snla datasets. The best fits among those 
parametrizations share the following common features 

• The inequality l|l.(jfl is violated for z ~ 0.3 ie the 
effective equation of state is w(z) < — 1 for z ~ 0.3. 

• The inequality 11.61) is respected for z ~ 0.3 ie the 
PDL is crossed at z ~ 0.2 — 0.3. 
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A representative easy to handle parametrization that 
shares these common features at best fit and has been 
one of the first introduced in the literature is the poly- 
nomial parametrization0 

q{z) = n 0m (l + zf + 02(1 + zf + ai(l + z) + a (3.1) 

where do = 1 — 02 — ai — i^om due to flatness. For a prior 
of fiom = 0.3 the parameter values that best fit the Gold 
dataset are a x = -4.16 ± 2.53, a 2 = 1.67 ±1.03. The 
effective equation of state parameter w(z) corresponding 
to the parametrization i|3.1|) is shown in Fig. 1 along with 
the corresponding form of H(z) at best fit. The form of 
H(z) corresponding to a cosmological constant with the 
same prior is also shown for comparison. In the redshift 
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FIG. 1: The effective equation of state parameter w(z) cor- 
responding to the parametrization 113. Ill along with the cor- 
responding form of H(z) — Hp2(z) at best fit. The form of 
H(z) = Hlcdm(z) corresponding to a cosmological constant 
with the same prior is also shown for comparison (dashed 
line). 

range < z < 0.25 where w(z) < —1 the inequality (|1.6|) 
is temporarily violated. The PDL is clearly crossed and 
therefore no minimally coupled scalar can account for the 
accelerating expansion crossing w(z) = —1 (see equation 
(|2.12l0 . We will thus examine if a non-minimally coupled 
scalar can account for this type of expansion. 

The parametrization i|3.1[l will be used for the recon- 
struction of the scalar tensor potentials along the lines de- 
scribed in the previous section. No attempt will be made 
to reconstruct the unique scalar tensor theory obtained 
from data. Instead we only examine if it is in principle 
possible to construct extended quintessence Lagrangians 
that lead to a crossing of the PDL in the way favored by 
the Gold dataset. We thus use arbitrary forms of $'(0) 
and first try the 'best case scenario' for the existence of 
F{z) setting &{z) = in (|2~TT}|) while using boundary 
conditions /(0) = 1, /'(0) = -2. The solution of (|2~TU|) 
for F{z) = (1 + z) 2 /(z) is easily found (see Fig. 2) to 
fulfil the requirement (|2.17H f i|2.18[l is imposed by the 



boundary conditions) for all redshifts z and is therefore 
observationally acceptable. 
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FIG. 2: The function F(z) can fulfil the requirements 12.171 - 
12.181 for a constant $(z) (eg $ = 0). 

We next ask the following question: 'Can the require- 
ments (I2.16f) - (I2.18II continue to hold if we use a small 
positive $'(z) 2 ?'. By numerically solving equation (|2.1U|) 
for a few trial functions <I>'(z) 2 and boundary condi- 
tions /(0) = 1, /'(0) = —2 it is easy to see that only 
small, decreasing with z functions <5>'(z) can accommo- 
date F(z) > for all redshifts z. Examples of such func- 
tions include 

$'(z) 2 = a (l + 2)" n (3.2) 
$'(z) 2 = ae- z (3.3) 

with a ~ 0.35 and n ~ 1. All such functions lead to 
qualitatively similar forms for f(z). 

In order to reconstruct the scalar potentials U(&), 
F(&) of the scalar tensor theory for the best fit q(z) of 
(|3.1|) we use the following steps 

1. Pick a function $'(z) = g(z) leading to F(z) > 
in the redshift range of interest. Such functions 
are given eg by equations l|3.2|l and l|3.3|l . Use this 
function and the best fit q(z) of equation l|3.1(l to 
solve equation (|2.1U|) numerically. 

2. Use the resulting /(z) along with the best fit q{z) 
in equation 12.9(1 to find U(z). 

3. To convert U(z) to U (<£•) solve the differential equa- 
tion <&'(z) = g(z) using the selected function g(z) 
with boundary condition $(0) = (set eg $ to at 
the present time). Use the solution $(z) to invert 
and find z($). 

4. Substitute z(<i>) in the U(z) found in step 2 and 
in F(z) = (1 + z) 2 /(z) found in step 1 to find 
U(z($)) = U($) and F(z($)) = F(<F). 
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Notice that the inversion z($>) can not be made for all 
values of $ but only for those values of $ appearing in 
the &(z) plot. This is a limitation because we have found 
that $(z) tends to a constant at large z (see Fig. 3) 
because of large cosmic friction at early times. 

The above steps have been implemented in a mathe- 
matica code available at 0] . For definiteness we have as- 
sumed a $'(z) of an exponentially decreasing form given 
in (|3.3() with a = 0.35. The resulting functions &(z), 
U(z) and F(z) are shown in Fig. 3 in the redshift range 
< z < 2. The potential U(z) is found to be a slowly 
increasing function of z while F(z) has a shallow mini- 
mum at z ~ 1 and increases monotonically beyond that 
redshift. The minimum of F(z) in Fig 3 (see also Fig. 
5) however, is deeper compared to that of Fig. 2 cor- 
responding to $'(z) = because F(z) has to balance, 
the attractive gravity of &'(z) ^ and therefore has to 
divert more from F(z) = 1. It is easily shown that the 
resulting potentials remain qualitatively similar for other 
acceptable forms of By inverting $(z) along the 
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FIG. 3: The reconstructed functions $(2:), U(z) and F(z) in 
the redshift range < 2 < 2. 

lines of step 3 the potentials U(Q) and F($>) may be 
found and are shown in Figs. 4 and 5 respectively. In 
Fig. 4 the numerically obtained potential corresponds to 
the thick dotted line. The existence of such potentials 
is our main result and not their particular form which 
is not unique given the form of H{z). Nevertheless for 
illustration purposes we attempt to fit the numerically 
obtained potentials with simple analytic functions. The 
continuous lines in Fig. 4 correspond to attempts to fit 
the numerically obtained potential with simple analytical 
functions. We have tried several simple analytical func- 
tions. A very good fit was provided by an exponential of 
the form 

U ($) = A e A * 2 (3.4) 

with A ~ 3. The fit was worse for other powers of <& in the 
exponential (the next best approximation U(&) ~ e A * 



(A ~ 5) is also shown in Fig. 4 and the fit is clearly not 
as good as with (|3.4|) '). 
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FIG. 4: The reconstructed potential t/(<3>) is well fit by e A * 
but not with other powers of $ in the exponent (eg e A *). 

It is straightforward to test our results, by substitut- 
ing the reconstructed functions U(z) and f(z) back in 
equation l|2.9() to solve for q(z) numerically (with bound- 
ary condition q(0) = 1) and confirm that we get back 
the best fit form l|3.1|l with the appropriate parameter 
values. This test revealed that the numerically obtained 
functions U and F of Figs 3-5 indeed lead to a prediction 
of exactly the best fit form of q(z). 
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FIG. 5: The reconstructed function F(&) has a local mini- 
mum at a small value of $ which is deeper than that of Fig. 
2. This minimum appears to be a generic feature found using 
other best fit parametrizations of H(z) as well. 

Alternative best fit forms of q(z) that violate the in- 
equality l| 1.6(1 crossing the PDL can be used to recon- 
struct potentials i ? ($), U($>). For example Ref. ^lj nas 
indicated a best fit q(z) parametrization that crosses the 
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PDL repeatedly, showing evidence for oscillations in red- 
shift space. This parametrization was shown to provide 
better fit to the Gold dataset than any other parametriza- 
tion tested in Ref. 11]. It is straightforward to recon- 
struct potentials F and U that reproduce the best fit 
oscillating H(z) parametrization and show that the ex- 
pansion history oscillations are inherited by the recon- 
structed potentials [43. Thus the U ~ e A * form of the 
reconstructed potential is not a particularly robust pre- 
diction of the Gold dataset. 



IV. CONCLUSION 

We have demonstrated that it is possible to reconstruct 
experimentally viable scalar-tensor potentials which pre- 
dict exactly the best fit parametrizations obtained from 
the recent Gold Snla dataset. These parametrizations 
share the common feature of temporarily violating the 
inequality H1.6(l or equivalently correspond to a dark en- 
ergy equation of state that crosses the PDL. This feature 
is impossible to reproduce in the context of single field 
quintessence or even phantom models. 

As an application we have used a quadratic polyno- 
mial parametrization for H{z) at its best fit with respect 
to the Gold dataset to reconstruct the scalar- tensor po- 
tentials F(<&) and J7($) taking into account consistency 
with solar system and Cavendish type experiments. The 
reconstructed scalar tensor potential U was found to be 
well fit by an analytical function of the form U($>) ~ e 3 * 
which may be thought of as physically motivated on the 
basis of SUGRA theoriesjl^lij- 

The derived potentials however are not uniquely deter- 
mined for two reasons 

• We used one observationally determined input 
function (H(z)) in the context of two coupled equa- 
tions to construct three output functions ( < &(z), 
F(z), U(z)). The additional solar system test ex- 
perimental constraints were not enough to lead to a 
unique determination of the three output functions. 

• The best fit form of H(z) depends on the type 
of parametrization considered. The polynomial 
parametrization considered for H(z) gives a very 
good quality of fit to the Gold dataset but is not 
unique. For example an oscillating parametrization 
gives a somewhat better fit, yet the reconstructed 
potentials are quantitatively different from those 



obtained from the polynomial parametrization at 
best fit. 

Thus our main result is that in contrast to minimally cou- 
pled quintessence, scalar tensor theories can reproduce 
the main features of the best fit Hubble expansion his- 
tory obtained from the Gold dataset. However, the pre- 
cise determination of the scalar tensor theory potentials 
requires more accurate Snla data and additional observa- 
tional input 152,153,0,11^113 which could come eg from 
weak lensing or large scale structure surveys providing 
the structure formation evolution history ^(z) = S m (z). 
Once S m (z) is known from observations its time evolution 
equation 



2H6 r , 



4irG cS {t) p 6 







(4.1) 



may be used in addition to the system (|2.9I) - (|2.10() to fully 
determine the scalar tensor Lagrangian. The unknown 
functions F(z),U(z) and $>(z) enter in equation l|4.1(l 
though the effective Newton's constant G e g(z) measured 
in Cavendish type experiments (Force = G c s m ^." 2 ) 
which is connected with the Newton's constant G en- 
tering in the scalar tensor Lagrangian (|2.1|1 by ^3] 



_ _ G (2ZF + A(dF/d<S>) 2 

G e S = = 



(4.2) 



F \2ZF + Z(dF/d<S>) 2 
In terms of the redshift z, equation (|4.1|) takes the form 



H 2 5" 



(q 



2V 



1 + z 



3, -Goff(£) 
2 ' — G — 



(4.3) 



The use of equation (|4.3|) to supplement the system (|2.9|) - 
(|2.10() and uniquely determine the functions F(z),U(z) 
and $(2) subject only to the observational uncertainties 
of the observed input functions H (z) and 8 m (z) is a par- 
ticularly interesting future prospect. 

The Mathematical!?} file including the numerical anal- 
ysis and the production of the figures of the paper can be 
downloaded along with the source files of the paper from 
the astro-ph archive or sent by e-mail upon request. 
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